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ABSTRACT 



Aims. We perform a deep search for galaxies in the redshift range 6.5 < z < 7.5, to measure the evolution of the number density of 

luminous galaxies in this redshift range and derive useful constraints on the evolution of their luminosity function. 

Methods. We present here the second half of an ESO Large Programme, which exploits the unique combination of area and sensitivity 

Q , provided in the near-IR by the camera Hawk-I at the VLT. We have obtained ~ 30 observing hours with Hawk-I in the K-band of two 

^ ■ high galactic latitude fields. We combined the K-band data with deep J and K Hawk-I observations, and with FORS1/FORS2 U, B, 

"^ ' V, R, I, and Z observations to select z-drop galaxies having Z - Y > 1, no optical detection and flat Y - J and Y - K colour terms. 

rrt ' Results. We detect 8 high-quality candidates in the magnitude range Y = 25.5 - 26.5 that we add to the z-drop candidates selected 

in two Hawk-I pointings over the GOODS-South field. We use this full sample of 15 objects found in ~ 161 arcmin^ of our survey 

to constrain the average physical properties and the evolution of the number density of z ~ 7 LBGs. A stacking analysis yields a 

best-fit SED with photometric redshift z = 6.85^q jj and an E(B - V) = 0.05 ^JJ^j. We compute a binned estimate of the z ~ 7 LF 

^ ' and explore the effects of photometric scatter and model uncertainties on the statistical constraints. After accounting for the expected 

^O , incompleteness through MonteCarlo simulations, we strengthen our previous finding that a Schechter luminosity function constant 

0^ ■ from z=6 to z=7 is ruled out at a >99% confidence level, even including the effects of cosmic variance. For galaxies brighter than 

CO ' Afi5oo = -19.0, we derive a luminosity density p^v = l-5^08 ^ lO^^erg i"' Hz' Mpc^^, implying a decrease by a factor 3.5 from 

VO ' z = 6 to z - 6.8. We find that, under standard assumptions, the emission rate of ionizing photons coming from UV bright galaxies is 

lower by at least a factor of two than the value required for reionization. Finally, we exploit deep Hawk-I / and K band observations 

i-^ ■ to derive an upper limit on the number density of M1500 S -22.0 LBGs at z ~ 8 (K-dropouts). 

^^ . Key words. Galaxies: distances and redshift - Galaxies: evolution - Galaxies: high redshift - Galaxies: luminosity function 
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1. Introduction ties r emains on the homogeneity (e.g. iMesinger & Furlanettol 

I2OO9I) and on the exact t imeline of the reionization process 



The search and study of galaxy populations at very high redshift Ve.g. iGaherani eTZI I2OO6I) . Whether the UV Ught emitted by 
IS one of the most promising research areas of today astrophysics gtai'-forming galaxies is capable of reionizing the Universe by 
and cosmology. It derives its importance on two different and ^j^^^^ ^^^^^^ ^^^^^^^ ^„ ^p^^ question that should be answered 
inten-elated aspects: 1) The estimate of the UV photon budget ^^^^^^^ ^^^ ^„^jy^i^ ^f j^j.gg ^^^^pj^^ ^f ^^^^ ^^^^^^^^ ^^-^^^^ 
provided by star-forming galaxies and its role on the reioniza- 
tion of the universe at z > 6; 2) The study of the formation and The search for high-redshift star forming galaxies has been 
the physical properties of the first bulding blocks of present-day carried out so far mainly with renditions of the Lyman Break, 
galaxies. or "drop-out" technique that has been proved to be extremely 



There is observational evidence that the Univers e is hig hly efficient at redshift fro i n 2 to 6 (e.g. i Steidel et al. 1995, 1999t 
ionized at z ~ 6 (e.g. iFan et all l2006t JTbtani et all I2006I) . in 'Ad elbereer etaP l2()04t iDicki nson et al. 2004; Giavalisco et aL 



agreement with the latest WMA P estimates of the Thomson op- 20041: lOuchi et al.1 l2004t iBouwens et al. 2007; McLure et aL 



tical depth (iKomatsu et alJ2010h . although significant uncertain- 120091) . or through naiTow-band studies targeting the Ly a emis- 

sion (e .g. live et al. 2006; Kashikawa et al. 2006; Ou chi et al.l 

Send offprint requests to: M. Castellano, e-mail: l2009bh . The application of the Lyman Break technique at z > 6 

castellanoOoa-roma. inaf.it has been performed, at first, in small areas with deep near-IR 
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J + H NICMOS data (e.g. iBouwens et aP |2004 . and it has 
recently acquired momentum thanks to the installation of the 
WFC3 camera onboard of the Hubble Space Telescope yield- 
ing to a sample of tens of faint Lym an Br eak galaxies (LBGs) 
(IBouwens etal.1 1201 Oct lOesch et al.l 120101: i McLure et al.ll20ld 
iBunker et alJ l2009t lYan et alJ l2009t IWilkins et al. 2010ab). 
In th e meantime, ground based surveys (Ouchietal. 2009a; 
IC^ak et al. 2009; Hickey et al. 2010; C astell ano et al. 2010, 
CIO hereafter ), along with refined an alysis of archival NICMOS 
observations dBouwens et alJl2010ah have expanded the number 
of bright LBGs known. 

The basic feature of the high redshift galaxy popula- 
tion that can be analysed through the present datasets is 
its UV luminosity function (LF). The current picture of the 
evolution of the UV LF points to a factor of 6-11 de- 
crease in the number density of UV bright galaxies from 
z ~ 3 to z ~ 6 (e .g. Stanwa v et al.l 120031 : IShimasaku et alJ 
120051: IBouwens et alj 12006.) . although some uncertainties are 
still present on the exact amount of evolut ion in the dif- 
ferent parameters of the Schechter function (Dickinson et al.' 
'2004; Giavalisco et al. 2004; Sawicki & Thompson 2006; 
Iwataetal. 2007; Yoshida etTaH l2006t IBouwens et alJ l2006l 
2007.: .Beckwithetal...2006) . At redshift above 6, most of the 



analysis indicates a strong evolution in the LF, mainly through a 
dimming of the characteristic magnitude M* and/or a decrease of 
the normalization factor (Bouwens et al. 2008; McLure et al . 
2010; Ouchi et al. 2009a; Yan et al. 2009; Castellano et alj2010t 
Bouwens et al. 20 10a) . The recent WFC3-based analysis by 
Oesch et al. (201(J) also found evidence for a steep faint-end 
{a ~ -1.8) , in agreement w i th the predictions of theoret- 
ical models (J Trenti et alJ 1 20101 : ISalvaterra et"an l 2010l ). LBGs 
searches around lensing clusters have also been performed, find- 
ing discrepant results that highlight the many challenges and 
uncertainties in thes e investigations (Richard et al. 2006, 2008; 
iBradlev et alj 120081: [Bouwens et all 120091: IZheng et all I2009D . 
Along with improved constraints on the LF at z ~ 7, the lat- 
est analysis of the WFC3 data have also provided a first estimate 
of the evolution at z ~ 8-9 that points to a further decrease 
in the LBG number density, and thus in the total amount of UV 
photons produced by young stars at these early epochs. 

The discrepancies among different works, both at z ~ 3 - 5 
and at z ~ 6 - 9 are most probably due to the effect o f cosmic 
variance (e.g. lTrenti & Stiavellill2008l:lRobertsonll201 0'). but also 
to the difficulties in avoiding systematic effects in the different 
estimates of completeness level, contamination from lower red- 
shift interlopers, vol ume elemen ts, and redshift distributions in 
the various samples dStanwav et al. 2008), all worsened by the 
known degeneracy among the parameters adopted to fit the LF. 

The strong decrease observed in the UV emission coming 
from relatively bright sources seems to imply that reionization 
cannot be explained on the basis of UV bright galaxies only. An 
increased number of low luminosity galaxies indicated by the 
steep faint end of the Schechter LF might play a decisive role 
in the reionization process. Large and reliable samples of high-z 
galaxies both at the bright and at the faint end of the LF are thus 
necessary to shed light on this issue, and, possibly to highlight 
the need to search for even more intriguing so urces of the reion- 
izing radiation wi th future facilities (see e.g. IVenkatesan et alJ 
12001 lMadau"erai...2004.) . 

Latest surveys have also given the opportunity of analysing 
the physical properties of high redshift galaxies, whose knowl- 
edge is also a decisive factor to understand the very role of 
these sources in the reionization process. Recent studies have 
given the first estimates of masses, ages and SFRs for sin- 




Fig. 1. Colour-composite image of the BDF (left) and NTTDF 
(right) fields, created using the weighted mean of Hawk-I Y, J 
and K images as red, the F0RS2 Z as green, and the weighted 
mean of FORS 1/FORS2 optical images as blue. 



gle z > 6.5 objects, the first constraints on the stellar mass 
density at these epochs, and have also raised an interesting 
debate on the possibilty that the first galaxes might be char- 
acterized by peculiar properties, like a very low dust con- 
tent, nearly primordial metallicity or top-heavy stellar initial 
mass fun ctions (Finkelsteinet al. 2009; Bouwens et al. 201()a; 
Gonzalez et al. 2010; Labb e et al.ll2010t ISalvaterra et all 1201 a 
.Schaerer & de Barros.2010.) . 

To give an answer to some of the a bove problems we 
are using the new VLT IR imager Hawk-I dPirard et al J 120041 : 
Casali et al. 2006; Kissler-Patig et al. 2008'), to conduct a deep, 
medium area survey in the Y band over four independent point- 
ings, aimed at the detection of relatively bright LBGs at 6.5 < 
z < 7.5. Thanks to the extreme efficiency and large field of 
view (7.5x7.5 arcmin) of Hawk-I, it is possible to easily reach 
Y ~ 26.5 AB at > 5cr (roughly corresponding to M1500 - -20.5 
at z = 7), over large areas in a reasonable amount of time (15 
hrs). 

In CIO we discussed the results of the first half of our sur- 
vey, covering a large fraction of the GOODS-S field, and we 
estimated a statistically significant (~99 % c.l.) decrease, with 
respect to z ~ 6, of the number density of UV-bright galaxies. 
In this paper we present the z ~ 7 candidates found in the sec- 
ond half of the survey, covering two other independent fields. 
We will constrain the evolution of the LF combining this new 
sample with the GOODS one. 

The paper is organised as follows. In Section 2 we present 
the imaging set and the multiwavelength catalogue; in Section 
3 the LBG selection criteria and the potential interlopers affect- 
ing the selection are discussed; in Section 4 we present our fi- 
nal sample of candidate z-drop LBGs. In Section 5 we discuss 
a stacking analysis of all the z-drop galaxies found in the four 
Hawk-I pointings, that are used to constrain the z > 6 UV LF in 
Section 6. In Section 7 we derive an upper limit on the number 
density of very bright z ~ 8 LBGs. A summary of our methods 
and results is provided in Section 8. 

Throughout the whole paper, observed and rest-frame mag- 
nitudes are in the AB system, and we adopt the A-CDM concor- 
dance model (//q = IQkm/s/Mpc, Q.m - 0.3, and Q,a = 0.7). 

2. Data 

2.1. Observations 

This work is based on deep F-band images obtained with the IR 
camera Hawk-I at the VLT, and on deep optical FORS2 obser- 
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Table 1. BDF - Observations 



Filter 


Instr. 


Exp. Time (s) 


Seeing (arcsec) 


Mag. Limit" 


V-High 


FORS2 




13800 




0.75 


29.1 


R-Special 


F0RS2 




11600 




0.63 


29.3 


I-Bessel 


FORS2 




4800 




0.70 


27.8 


Z-Gunn 


FORS2 




64800 




0.59 


28.6 


Y-Open 


HAWK-I 




56940 




0.52 


28.3' 


J-Open 


HAWK-I 




18720 




0.54 


26.5 


Ks-Open 


HAWK-I 




30060 




0.44 


26.0 






a 


-S/N=l 












b 


- Y=26.5 at S/N= 


=5 





Table 2. NTTDF - Observations 



Filter 


Instr. 


Exp. Time (s) 


Seeing (arcsec) 


Mag. Limit" 


U-Bessel 


FORSl 




32876 




0.84 


27.8 


B-Bessel 


FORSl 




16064 




0.56 


28.9 


V-Bessel 


FORSl 




10500 




0.47 


29.0 


R-Special 


F0RS2 




14000 




0.79 


28.4 


I-Bessel 


FORS2 




7830 




0.61 


28.0 


Z-Gunn 


FORS2 




46386 




0.60 


28.4 


Y-Open 


HAWK-I 




54180 




0.49 


28.3* 


J-Open 


HAWK-I 




14400 




0.47 


26.7 


Ks-Open 


HAWK-I 




24720 




0.39 


26.3 






a 


-S/N=l 












b 


- Y=26.5 at S/N= 


=5 





vations. We use data collected through a dedicated ESO Large 
Programme in 2008 and 2009. The first set of data, covering 
two adjacent regions of the GOODS-S field has been presented 
in CIO. Here we present the analysis of two other pointings 
(Fig.[T]i, chosen for the wealth of deep, public observations pre- 
viously exploited by other authors to search f or z ~ 4 - 6 LBGs: 
the B DF field at Ra=336.98°, Dec^-35.17°(lLehnert & Bremeij 
l2003l) . and the New Technol ogy Telescope Deep Field ( NTTDF) 
at Ra=181.36°, Dec^-7.72° (lArnouts et alj[T999t iFontana et alj 
[2000. 2003) . The total exposure time is 15h49m for BDF and 
15h03m for the NTTDF in the Y band. 

The Y band images were reduced using standard techniques 
for IR data - flat fielding, sky subtraction among consecutive 
frames, and final coaddition. The reduction procedure, which is 
described in detail in our first paper CIO, has been specifically 
designed to enhance the reliability of the images at the faintest 
fluxes, and to get rid of persistence effects and cross-talk reso- 
nances. 

We determine an FWHM of 0.52 +0.01 arcsec (^ 4.9 pix- 
els) in the final coadded BDF image and 0.49 ±0.01 arcsec (^ 
4.6 pixels) in the NTTDF one. Image zeropoints were computed 
using the standard stars observed during the same night and at 
similar airmasses. Reference fluxes were converted to the pho- 
tometric system and filter set used in this paper, as described in 
CIO. 

We obtained the absolute rm.s. maps for each pointing 
by computing the r.m.s. in each individual image (using the 
Poisson statistics and the instrumental gain) and propagating 
self-consistently this rm.s. over the whole data reduction pro- 
cess. The typical 5cr magnitude in one arcsec^ is in the range 
26.7-26.8 over more than 60% of the whole image, and > 26.2 
in 85% of the image - the rest of the images being shallower 
because of the gaps between the four Hawk-I chips. 

A wide wavelength coverage is needed to reliably select high 
redshift LBGs excluding lower redshift interlopers and red and 
dusty galaxies at intermediate redshift. To this aim we obtained 



or re-reduced deep observations of both fields ranging from the 
blue to the near-IR, matching all images to the Y band pixel-size 
and astrometric solution. 

Along with the main F-band pointings, we also acquired 
deep J and Ks Hawk-I observations of both fields. We also 
obtained ~ 7 hours of FORS2 Z band coverage for each 
field, that we coadded with the already existing FORS2 images 
(Fontanaet al. 2003) to reach the required depth. We also re- 
reduced the archive U, B, V, R, I FORS2 and FORSl observa- 
tions of the NTTDF, and the FORS2 R and / images of the BDF 
Finally, we obtained ~ 4 hours of y-FORS2 observations on the 
BDF. The full dataset is presented in Tab. [1] and Tab.|2l 

2.2. The photometric catalogue 

2.2.1. Detection 

We obtaine d the photometric catal ogue using the SExtractor 
code V2.5 (iBertin & ArnoutsI [T996h and the Y band as detec- 
tion image with the r.m.s. map derived as described above. Since 
high redshift galaxies are almost unresolved in ground-based im- 
ages, and SExtractor's MAG_BEST are known to underestimate 
the total flux of faint objects (Y > 24 in our case), we chose 
to use aperture-corrected total magnitudes. We computed aper- 
ture magnitudes in a 2 FWHM diameter and coiTected them to 
total magnitudes adopting aperture corrections from bright non- 
saturated stars in each field. While this choice might give slightly 
underestimated fluxes for the more extended high redshift can- 
didates, we can easily take into account this systematic through 
the simulations that we use to estimate the LF (Sect. 16.1b that 
are based on the observed profile of LBGs with know n spec- 
troscopic redshifts 5.5 < z < 6.2 (IVanzella et al.l 120091) in the 
GOODS-S ACS images. 

We optimised the SExtractor parameters involved in the de- 
tection process through the analysis of a 'negative' image as dis- 
cussed in CIO, adopting the set of parameters that minimises 
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the ratio between 'negative' and 'positive' detections at the faint 
end of the number counts. As expected, we find that the best 
parameters for faint objects detection on BDF and NTTDF Y- 
band images are the same adopted for the similar set of images 
over GOODS-South: we require 10 contiguous pixels each at 
S /N > 0.727, coiTesponding to a 2.3cr detection, and we re- 
strict the analysis to the regions where the r.m.s is less than ~ 
1.5 times the lowest value. With this choice of parameters, we do 
not find any detection on the negative images at F < 26.2, and 
a fraction of negative detections less than 5% of the real ones 
at fainter magnitudes. However, a posteriori, the latter value 
overestimates the actual rate of spurious detections. Indeed, all 
spurious sources should appear as "drop-out" candidates with a 
single-band detection. On the contrary, all the Y > 26.2 objects 
in our z-drop sample are confirmed by detections in other IR 
bands. Indeed, as we discuss also in CIO, the test on the negative 
image is probably influenced by non-trivial issues concerning 
the subtraction of the background or a potential asymmetry in 
the noise distribution. 



2.2.2. The multicolour catalogue 

A multiwavelength catalogue containing self-consistent magni- 
tudes in all available bands was built running SExtractor in dual 
mode using the F-band Hawk-I image as the detection image 
with the detection parameters indicated above. Aperture fluxes 
were computed within a 2FWHM aperture and converted to to- 
tal applying appropriate aperture corrections in each band. 

The typical Icr limiting magnitudes in a 2FWHM aperture 
are in the range 27.8-29.3 for the optical bands, J ~ 26.5-26.7, 
andKs ~ 26.0-26.3. The coiTesponding Icr limiting magnitude 
in the Z band, which is used to define the 'dropout' selection, is 
~ 28.6 in BDF and ~ 28.4 in the NTTDF (see Tab. [Hand Tab. 
Ell. 

For each field we defined the total areas where the image 
depth is sufficiently homogeneous. The candidates found in this 
area will be used for the evaluation of the LF We used F-band 
detected objects only in the regions selected on the basis of the 
negative image test explained above. In addition, we also masked 
borders, CCD defects and noisiest regions in the other images of 
our data-set. The areas selected in this way correspond to ~ 7 1 % 
of the y-band coverage in the BDF, and ~ 56% in the NTTDF 
(due to strong vignetting in the Z-band image). As a result, the 
total area used for z-drop detection amounts to 71.7 arcmin^. 
We will subtract to this value the fraction of area covered by 
lower redshift objects (~9%) to estimate effective volumes in 
Sect.|6]andSect.|71 



3. The selection of z> 6.5 galaxies 

3. 1 . The dropout criterion 

We select candidate z > 6.5 galaxies using the "drop-out" tech- 
nique adapted to our filter set and imaging depth. 

In order to individuate the appropriate selection criteria we 
estimated the expected colours of high-redshift star forming 
galaxies (black points in the right panel of Fig. |2]l on the basis 
of the models of Chariot and Bruzual 2007 (.Bruzual.2007aabh 
with the same range of free parameters as in CIO: MetalUcity: 
0.02, 0.2 and 1 Z©; age from 0.01 Gyr to the maximal age of the 
Universe at a given z; E(B-V)=0...0.2 (Calzetti et al. 2000). Ly-a 
rest-frame equivalent width in the r ange 0-200 A. Intergalactic 
absorption following iMadatJ (Il995h . The same range of model 



sooo 




6 6.5 7 

Hedshifl 



Fig. 2. Left: Z-Y colour of star forming galaxies as a function of 
redshift. In the upper part, the efficiency curve of the two filters 
is shown, computed at observed wavelength of a Lyman-c emis- 
sion at the corresponding redshift. Right: Z - Y vs. Y - J colour 
diagram showing the expected colours of LBGs (same as in left 
panel, black points), passively evolving galaxies (red squares) 
and reddened starbursts (green cir cles) at 1 .5 < z < 4 and cool 
dwarf stars from the templates of iTsuii et al.l (l2004t) (magenta 
stars). Galaxy colours are computed according to CB07 models, 
see text for details on the adopted parameters. In both panels 
lines indicate the relevant colour selection criteria discussed in 
Sect. [3] 



parameters will be used as baseline for the MonteCarlo simula- 
tions used to estimate the LF in a self consistent way, see Sect. 

m 

As shown in the left panel of Fig. |2] galaxies at z > 6.5 show 
an increasing Z-Y colour which is due to the sampling within 
these two filters of the sharp drop shortward of the Lyman-a, 
where most of the photons are absorbed by the intervening HI 
in the intergalactic medium. The drop in the flux observed short- 
ward of the Y band is analogous to the one used to select star 
forming galaxies at lower redshifts, like /-drops at z ~ 5, V-drops 
at z ~ 4 etc.: the major difference with respect to the standard 
Lyman break technique being that the Y band does not sample 
the continuum around 1500 A but a region shortward of it, con- 
taminated by both the larger IGM absorption at z > 6 and by 
the Lyman-a emission line. These eff'ects can only be accurately 
accounted for by realistic imaging simulations, as we discuss in 
detail in CIO and in section 16. II of this paper. Following this 
test, we choose Z - Y > 1 as our main selection criterion to 
select z > 6.5 galaxies. Given that the Z-band observations, as 
well as the optical ones, used in the present paper are slightly 
shallower than the GOODS-ACS ones, we limit our selection at 
Y < 26.5 instead of the Y < 26.8 adopted in CIO. 

The selection of z-drop galaxies cannot be solely based on 
the Z-Y colour, since other classes of objects can display a red 
Z-Y colour similar to that of z > 6.5 galaxies. Selection criteria, 
both in the optical and in the IR bands, are thus necessary to 
individuate a reliable sample of z ~ 7 galaxies. 

3.2. IR colour selection 

We tailored our IR colour selection to exclude any possible con- 
tamination in our z-drop sample from known classes of lower 
redshift objects: 

i) We modelled passively evolving galaxies and dusty star- 
burst galaxies at z > 1.5 with a suitable set of spectral syn- 
thesis models. We use the same CB07 library as for high red- 
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shift galaxies to predict the colours of such objects at 1.5 < 
z < 4, using a combination of short star formation exponential 
timescales (0.1 - 1 Gyrs) and ages > 1 Gyr to reproduce pas- 
sively evolving galaxies, and con stant star-form ing models with 
0.5 < E(B-V) < 1.5 (adopting a FCalzetti et aill200 0. extinction 
law) for the dusty starbursts. As shown in Fig. |2] (right panel) 
these galaxies also show a large IR colour term. 

To exclude these objects, we adopt the same additional crite- 
ria on IR colours as in CIO : (Z-Y)>(Y-K); (Z-Y) > 0.5-i-2.0(Y- 
J); (Y-J)<1.5;(Y-K)<2.0. 

ii) Cool (Teff < 1500 K), low-mass stars, and substellar ob- 
jects of the T spectral class have infrared spectra that are dom- 
inated by the CH4 a nd H2O absorption bands and by Ho reso- 
nant absorption (e.g. lChabrier et al.ll2005l : lBurgasser et al.ll2006l) 
that produce a sharp break in their IR colours. We used the 
most up-to-date estimate of the T-dwarfs number density (as ob- 
served in the J band) of Burgasser et al 2007 to compute the 
expected number of faint, cool dwarfs in our fields. Adopting 
an average Y - J colour of 0.8 mags estimated from the cata- 
logue of observed dwarfs compiled by'Leggett et al.' (201(f), and 
consid ering the dependence on galactic latitude as in Burgasser 
(l2004l) . we estimated that ~ 0.6 stars of spectral types T0-T8 
with Y < 26.5 are expected to fall in each of our fields. However 
the exact number of expected contaminating dwarfs depends on 
the still uncertain parameters constraining the IMF and their spa- 
tial distribution inside the disk and the halo of the Galaxy: a 
pessisimistic estimate (see CIO) gives a nearly double surface 
density of cool dwarfs in our point ings. For th i s reaso n we used 
the synthet i c spec tral libraries by iTsuii et al] (|2004|) (see also 
lTsuiill2002l l2005h to check whether it is possible to define se- 
lection criteria discriminating between LBGs and cool dwarfs. 
As shown in the right panel of Fig. |2] cool dwarfs appear redder 
than z-dropouts in the Y - J colour, and the Y - J criterion we 
adopt allows us to exclude these objects from our selection win- 
dow. We note that the brown dwarf discovered in the NTTDF by 
iCubv et al. (1999) . having Z -Y ~ 2 and no optical detection, is 
consistently excluded from our high-z sample on the basis of its 
Y - J colour. 

iii) Finally, we cross-checked each object selected according 
to the above criteria against variability, by analysing images ac- 
quired at different epochs. The BDF observations have been split 
in two separated epochs with a 3 months gap (September and 
December 2009), while the NTTDF have been observed during 
four runs in January, February, April and May 2009. We veri- 
fied that all the objects in our sample are clearly detected, and 
that they have a consistent total flux (within 2cr), in the different 
epochs. In the NTTDF case, we checked that a detection > 3cr of 
the faintest candidate {Y ~ 26.5) was possible in the two epochs 
with larger integration time. 

We summarise here the full set of colour selection criteria: 

Y < 26.5 

(Z-Y) > 1.0 

(Z-Y) > (Y- K) 

(Z-Y) > 0.5 + 2.0(Y - J) 

(Y-J) < 1.5 

(Y-K) < 2.0 



Team, in preparation) to select z-drop galaxies and to exclude 
lower redshift interlopers showing significant detection in the 
optical bands. The main concern we have to consider to pro- 
vide a z-drop selection as clean as the one in the GOODS field 
regards the difference in resolution between FORS2 optical ob- 
servations of BDF and NTTDF and their corresponding ACS- 
GOODS images we used to remove interlopers from the colour- 
selected sample. 

Indeed, in CIO we found that a sample of galaxies selected 
with IR criteria only is populated also by faint contaminants 
showing significant detection in filters covering wavelengths 
shorter than the redshifted Lyman limit at z > 6 (U, B, V, R, I) 
where high redshift LBGs are not expected to present any flux. 
These objects are, in most cases, clearly extended, but their spec- 
tral energy distributions cannot be reproduced by a straightfor- 
ward application of the CB07 models. While determining their 
nature is beyond the scope of the present analysis, we note that 
they might be faint galaxies with a very blue continuum whose 
SED is altered by strong emission lines such as in unobscured 
AGNs, or in star-forming galaxies like the blue compact dwarf 
galaxies (Izotov et al. 2004, 2007) or the ultra strong emission 
line galaxies (USELs, Hu et al. 2009). Potential contamination 
of z~ 7 samples due to an unknown class of objects with no op- 
tical d etection (< 2cr) has also been suggested bv ICapak et al.1 
(2009). Their objects are brighter than those found in our fields 
but display similar colours. Given the unknown nature of this 
contaminants, at present, the only feasible approach is to adopt 
more stringent criteria on the optical non-detections. Follow- 
up spectroscopy of z~ 7 candidates and faint contaminants is 
needed to fully evaluate the impact of this population on high 
redshift studies. 

In our analysis of the GOODS-South field we adopted very 
strict selection criteria in the blue bands in order to exclude these 
contaminants, measuring the S/N ratios in small apertures (0.6") 
exploiting the high resolution of ACS images. 

In order to obtain optical selection criteria as effective as 
the ones used with the GOODS dataset, we performed tests 
computing S/N ratios and photometry on the GOODS-ACS im- 
ages degraded and smoothed to the depth/seeing of the BDF 
and NTTDF corresponding ones. We then re-selected GOODS 
dropouts on "mock" BDF/GOODS and NTTDF/GOODS cat- 
alogues built in the same way as the real BDF and NTTDF 
catalogues. We verified that the criteria already adopted in the 
GOODS fields are effective in the NTTDF case: S /N < 2o-s/n 
in all the optical bands and < Ios/n in at least four of them. 
In the BDF, given the absence of U and B images and the 
slightly shallower / imaging, we adopted the conservative cri- 
terion S /N < IcTs/N in all the optical bands. We verified that 
this criterion allows us to safely remove all those objects, up 
to y=26.5, that have been verified to be lower redshift contam- 
inants on the basis of GOODS-ACS and, whenever possible, 
UDF-ACS photometry. 

The parameter cts/n indicated above is the r.m.s. of the S/N 
distribution estimated, as in CIO, dropping random apertures in 
portions of the images free of detected objects. This procedure 
allows us to take into account the sky noise distribution and 
the presence of faint, undetected foreground objects at the same 
time. 



3.3. Comparison with fhe GOODS-ACS dataset 

In our analysis of the GOODS-South field we exploited the ACS 
V2.0 B, V, I, Z observations (M. GiavaUsco and the GOODS 



4. Detected z> 6.5 galaxies 

Adopting the selection criteria outlined above we find a total of 
eight candidates, three in the BDF and five in the NTTDF field, 
whose coordinates, Y magnitudes and Z-Y colours are listed 
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Fig. 3. Thumbnails showing the images of the 8 selected high-redshift candidates in the different observed bands. 



in Tab [3] Thumbnails of the candidates are presented in Fig. |3] 
We note that two of them are clearly detected, and two others 
are marginally detected (S/N~ 2), in the Z band. Three of the 
five candidates present in the NTTDF are also detected at S/N~ 
2 - 4 in the J and K bands, thanks to the slightly deeper images 
available for this field (see Tab |2]i. We also verified that each 
candidate is undetected in the image obtained as the weighted 
sum of its V, R and / observations. 

As a final check we performed a stacking of all the objects 
in the available images. This test allows us to confirm the non- 
detection in the optical images, and to obtain a clear detections 



in the J (S/N~ 5) and K band (S/N~ 4) stacked images. The 
stacked object shows an average colour Z - Y ^ 1.6. 

In the following sections we will combine this sample z- 
drop candidates found in the BDF and NTTDF pointings with 
the sample discussed in CIO, obtained from the two point- 
ings over the GOODS-South field, to find their average prop- 
erties through a stacking analysis, and to constrain their LF The 
GOODS z-drop sample includes seven candidates in the range 
Y ~ 25.5-26.7, selected through colour selection criteria analo- 
gous to the ones outlined above, whose reliability has also been 
checked on the available IRAC and NICMOS observations. 
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Fig. 5. Best-fit SED to the stacked photometry, with relevant photometric redshift at z = 6.85. 
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Fig. 4. Thumbnails showing the stacked images of the 15 high- 
redshift candidates selected in the GOODS 1, GOODS2, BDF 
and NTTDF fields. The observed bands are shown in the legends. 

Table 3. Candidates in the BDF and NTTDF fields 



ID 



R.A. (deg) DEC. (deg) 



Z-Y S/N (Y) 



BDF^21 
BDF_3299 
BDF_5905 
NTTDF.1479 
NTTDF.1632 
NTTDF.1917 
NTTDF.6345 
NTTDF.6543 



336.9444 
337.0511 
337.0230 
181.3429 
181.3857 
181.3212 
181.4039 
181.3834 



-35.1188 

-35.1665 

-35.2094 

-7.6813 

-7.6835 

-7.6877 

-7.7561 

-7.7595 



25.86 
26.15 
26.24 
26.12 

26.47 
26.32 
25.46 
25.75 



2.13 
>2.4 
1.20 
1.97 
1.61 
1.58 
1.45 
>2.6 



10.2 

7.8 

7.6 

8.4 

6.1 

7.1 

15.6 

12.0 



didates in each image. The stacked object shows an S /N > 5 
detection in the Z, J and K bands and a non-detection in all the 
optical bands, corresponding to an {optical - Y) colour of > 4 
magnitudes. We use the magnitudes estimated for the stacked 
object to find the photometr ic redshift and physical parameters 
through our photo-z code ( Giallongo et alJfl998l : iFontana et al.l 
[2000) exploiting a x^ minimisation procedure to find the best- 
fitting spectral template to the observed colours among the full 
CB07 library. While the ACS optical filters used in the GOODS 
field have different passbands with respect to the FORS2 ones 
used in BDF and NTTDF fields, this is not a significant con- 
cern since they all span a wavelength range where no flux is 
expected for z > 6 objects. In turn, the small difference be- 
tween FORS2 and ACS Z-band filters does not provide signif- 
icant variations in the redshift selection window defined by the 
Z-Y colour which is the main constraint to the photometric 
redshift. The resulting SED provides a unique photometric red- 
shift solution at z = 6. 85;^yj*J. Relevant thumbnails and SED 
are shown in Figure |4] and \5\ Given the absence of IRAC, 
most physical parameters are largely unconstrained, apart from 
the E(B-V) parameter whose estimate is mostly based on the 
_Y - J and Y - K colours. We find that our stacked SED is fit- 
led by an E{B - V) = 0.05+°^^ at a 68% confidence level. This 
value is consistent with the E(B -V) distribution obtaine d from 
the analysis of z~ 7-8 objec ts bv lFinkelstein et al] (l2009t) and by 
ISchaerer & de BarrosI (12010'). Our best fit E(B-V) indicates a low 
dust content for z ~ 7 galaxi es, in agreemen t also with the best- 
fit Ay values found by Gonz alez et al.l (1201 Oh and bv lLabbe et al.l 
( 1201 Oh for the mean SED of their z-dr op samples, an d with the 
""blue UV continuum slope measured bv lBouwens et al.. (■2010d') . 



5. Mean properties of Hawk-I z ~ 7 galaxies 

We perform a weighted mean of the images in the available fil- 
ters from the U to the K band for all the 15 objects detected 
in the Hawk-I fields. We did not attempt a similar stacking of 
the IRAC images, since most of our GOODS candidates are par- 
tially or extremely blended with other foreground sources, and 
the candidates in the other fields are either not covered by IRAC 
observations, or they are present in shallower exposures with re- 
spect to GOODS. We matched the ACS images to the Hawk-I 
PSF and masked all the foreground objects surrounding the can- 



6. Tlie evolution of the LF 

6.1. MonteCarlo simulations 

When small galaxy samples are used to constrain the high- 
redshift LF, it is necessary to exploit detailed imaging simula- 
tions to appropriately treat the systematic effects arising from 
faint object detection, and from the application of colour selec- 
tion criteria. To this aim we use the CB07 synthetic libraries de- 
scribed in Sect. [3] to produce, for each field, a set of ~ 8 x 10^ 
simulated LBGs with redshift in the range 5.5 < z < 8.0 and 
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Table 4. Stepwise determination of the UV LF 



Mag. Range 



<p (10-'^ M pc-^ mag-^) 



-21.4<Mi5oo 
-20.8 < M,5oo 



< -20.8 

< -20.0 



0.39: 
1.81: 



0.20 
0.54 



observed magnitudes computed in the same filter set used for 
the observations. These galaxies are placed at random positions 
of the F-band images, and catalogs are extracted exactly as in 
the original frames. To avoid an excessive crowding in the simu- 
lated images, we include only 200 objects each time, after mask- 
ing the regions of the images where real objects are present. As 
in CIO, we randomly assign to each of our simulated galaxies 
the light profile of one of the four most distant spectroscopically 
confirmed LBGs ob served with ACS in GOODS (z = 5.5 - 6.2, 
IVanzella et al.l2009l) . after convolving it with the relevant Hawk- 
I PSFs. 



6.2. Stepwise LF 

The magnitude range covered by our survey, Y ^ 25.5 - 26.7 
roughly coiTesponds to the UV continuum magnitude range at 
M < Mt. For this reason, we first perform a binned estimate 
of the number density of the Hawk-I z-drop galax ies through 
the stepwise method (see, e.g. Bouw ens et al.ll2008!) . The step- 
wise estimate is a non-parametric method based on the assump- 
tion that the rest-frame LF of galaxies can be approximated by a 
binned distribution, where the number density 0, in each bin is a 
free parameter. To evalute also the potential systematics and the 
effects of observational uncertainties in this kind of estimates, 
we use two different procedures to compute the stepwise LF. 
The first one is the procedure commonly adopted in the Utera- 
ture based on the average relation between the observed Y and 
the UV continuum magnitude at 1500A (Misoo), and on an es- 
timate of the completeness in the different UV magnitude bins. 
The second, more conservative, procedure takes in consideration 
the uncertainties in the Y-Mism conversion due to photometric 
scatter, to the redshift distribution and to the intrinsic proper- 
ties of different galaxy models. In a separate work we will com- 
bine this stepwise analysis with similar estimates at fainter and 
brighter magnitudes to determine the Schechter parameters at 
z ~ 7 in a self-consistent way (Grazian et al. 2010, in prepara- 
tion). 



6.2.1. Stepwise LF from the average Y-UV relation 

Through a linear regression we compute the average Y - M1500 
relation at the median redshift of our sample (z=6.8) for the 
CB07 models of z-drop galaxies. We then divide our sample in 
two bins centered at M1500 - -2 LI and Misoo = -20.4, and 
use the imaging simulations to estimate the completeness of our 
selection. Finally, we convert the redshift dependent complete- 
ness distribution into effective volumes of our survey at these 
magnitudes. The values of the stepwise LF estimated in this way 
are reported in Tab. |4] and plotted as filled squares in Fig|6] with 
vertical error bars given by Poisson uncertainties in the number 
counts. The horizontal error bars indicate the relevant magnitude 
range of each bin. 




Fig. 6. Number densities in two rest-frame magnitude intervals 
estimated for our Hawk-I data set in a stepwise form with a stan- 
dard F-UV conversion of the observed number counts as dis- 
cussed in the text (black filled squares), or with a x^ method 
considering also photometric and model uncertainties (black 
empty circles). Other points a re from [Bo uwens et aP (1201 Oah 
(NICMOS, red empty squares), 'Ouchi etal. (20 09ai) (SUBARU 
blue empty squares and upper limits) and Oe sch et al] (1201 Oh 
(WFC3-UDF, magenta empty circles). For a comparison we 
sho w the recent determinations of the LF at z ~ 6 
by Bouwens etal. (l2007l) (B07, green dot-dashed line) and 
iMcLure et al.1 ('2009') (M09, blue dashed line). The black solid 
line is the best-fit LF obtained by combining the stepwise points 
shown in the figure with new determinations of the binned den- 
sities from WFC3-UDF and WFC3-ERS data (Grazian et al. in 
preparation). 



6.2.2. Introducing photometric uncertainties in the Stepwise 
LF 

A more conservative estimate can be computed assuming a 
stepwise LF made of three bins in the wider magnitude range 
-22.0 < Misoo < -19.0. This interval takes into account the 
photometric scatter and the variation of the Y - Misoo relation 
with redshift and galaxy models (see Fig|7]i. We assume a fixed, 
constant, reference density 0,,./, and we exploit the set of simu- 
lations described in Sect. 16.1 I to compute for each field the dis- 
tribution of observed magnitudes originating in each rest-frame 
bin for LBGs in the redshift range sampled by our colour selec- 
tion. The simulated number counts are then scaled to the rele- 
vant observed areas and summed together. Finally, we find the 
combination of binned densities 0,- - w, ■ (pref, that best repro- 
duces the total number counts of our survey, where w, are mul- 
tiplicative factors to the reference density that we determine by 
comparing observed and simulated distributions through a sim- 
ple ;^'^ test. We plot as black empty circles in Fig|6]the two bins at 
Mi5oo < -19.8. The third, faintest, bin at M1500 > -19.8 yields 
only a conservative upper limit and it is not represented in the 
figure, but it is anyway necessary in this procedure to consider 
the effect of Malmquist bias. Vertical error bars indicate the sta- 
tistical uncertainties given by the x^ test. 
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Fig. 7. The normalized distribution of UV continuum magni- 
tudes (estimated from the average Y - M1500 relation) for the 
15 Hawk-I candidates divided in two bins (dashed histograms). 
The solid curves show the expected distributions of UV magni- 
tudes for objects in the same observed ranges when photometric 
uncertainties are taken into account through MonteCarlo simu- 
lations. 



The two methods give consistent results, and they are in per- 
fect agreement with other stepwise estimates in the same mag- 
nitude range (see Fig. |6]l. However, the error bars and the rele- 
vant magnitude range are much larger when using the x^ mini- 
mization procedure. While an average conversion from observed 
to rest-frame magnitudes, along with an estimate of effective 
volumes, can provide a first order-of-magnitude estimate of the 
binned number density of LBGs, we emphasize that significant 
statistical uncertainties can arise due to photometric scatter, and 
to the different relation between Y and UV continuum magni- 
tudes for different galaxy models and redshifts. 

6.3. Maximum Likelihood LF 

We estimate how significant is the evolution of the LF at z > 6 
adopting a maximum likelihood approach. This method allows 
us to compare the observed number counts to those predicted 
for different evolving Schechter LFs (Schechter"l97d) after ac- 
coun ting for the expected systematics in the detection process 
(e.g. [Bouwens et al.ll2007l: iMannucci et alj|2007l: iMcLure et alj 



l2009l) . As in CIO, we assume that the LF can be described by 
a Schechter function with parameters an d M, evolving from 
their value at zo =6.0 (iMcLure et al.ll2009h according to the fol- 
lowing parametrisation: 

logi^iz)) = logi4>{zo)) + dlog{(l>)/dz ■ iz - zo) 

M,(z)^M,(zq) + m:-(z-zo) 

Since our faint limit is close to the expected value of the char- 
acteristic luminosity M,, we fix the faint end slope to the value 
a = -1.71 of the z ~ 6 LF bv'McLur e et al.l (l2009l) . We explic- 
itly tested that no appreciable differences are found when fixing 
a to different values (a - -1.4, -2.0). For a broad range of val- 
ues of the evolving terms M'^ and dlog{(p)/dz (see Fig.[8]l we sim- 



X^ / n.d.f. 
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Fig. 8. ;t'^ contour levels for the dlog(^)/dz. Ml parameters de- 
rived for the Schechter-like LF considering all the four Hawk-I 
fields. The lower and left axis refer to the evolutionary terms 
M ' and dlog((t>)/dz wit h respect to the best-fit z=6 parameters 
of iMcLure et al.l ( l2009l) (black point and errorbars). The upper 
and right axis refer to the M, and <p values at the median red- 
shift estimated for our sample (z=6.8). Grey points a nd errorbars 
mark th e positi on of the z ~ 7 b est fit parameters bv lOuchi et al.l 
(I2009al) and by iBouwens et alj (|2008). The black solid fine in- 
dicates the 99% c.l. region estimated on the basis of the two 
GOODS-South pointings only (CIO). 



ulate, under a MonteCarlo approach, the redshift z and UV mag- 
nitude Mi5()o for a population of 3 x 10^ galaxies. These objects 
are randomly extracted from the larger database of simulated 
galaxies described in Sect. 16.11 which encompasses a broad 
range of the physical parameters determining the rest frame pho- 
tometry, like E(B-V), metallicity, Lya EW etc. 

The distributions of Y magnitudes and Z- F colours for each 
simulated population are scaled to the observed area in each of 
the fields and compared to the observed ones with a maximum 
likelihood test, under the assumption of simple Poissonian statis- 
tics. For each of the two distributions, and for each field, we build 
the likelihood function X: 



i-U' 






(1) 



where Nobsj is the observed number of sources in the magnitude 
(colour) interval /, Nexp,i is the expected number of sources in the 
same magnitude (colour) interval, and 11, is the product symbol. 
For each field, we associate to every model a likelihood com- 
puted as the product of those obtained for the magnitude and 
colour distributions separately. We then compute a final likeli- 
hood as the product of the GOODS, BDF and NTTDF likeli- 
hoods. 

The colour plot in Fig. [8] shows the 68%, 95%, and 99% 
likelihood intervals on the evolutionary terms Ml and dlog{(p)/dz 
(left and bottom axes) and for the resulting Schechter parameters 
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at the median redshift z=6.8 of our sample (top and right axes) 
for the combination of all the four Hawk-1 fields. In the same 
plot, the colour code refers to the x^ distributio n obtained u nder 
the usual assumption x^ - -2.0 ■ /n(X) (e.g. ICashl[l979h . We 
reject at > 99% confidence level the hypothesis that the LF re- 
mains constant in both parameters above z—6 (dlog((fi)/dz — 
and M[ - dM^/dz - 0, black point in Fig. HI). In Fig. [8] it 
is also shown the 99% c.l. region on the Schechter evolutionary 
terms estimated on the basis of the two GOODS-South pointings 
only (CIO). Although the degeneracy between M» and (p is still 
present, the analysis of the BDF and NTTDF fields considerably 
reduces the allowed parameter space. 

The region of allowed values for the LF parameters in our 
final likelihood map points to a pronounced decrease of (f> along 
with a brightening of M, with redshift. However, the best-fi t val- 
ue s for M, and (p at z ~ 7 derived by 'Ouchi et al.' ('2009a') and 
bv lBouwens et al.l (l2008h (see also|Bpuwens et al. 2010b), indi- 
cating a constant or slightly dimming M», still fall within the 2cr 
region constrained by our maximum likelihood (grey points in 
Fig. [8]), and they are consistent with our estimate once the un- 
certainties are considered. We argue that cosmic variance (see 
Sect. 16.41 ) and the limited sample of very bright objects available 
may explain the discrepancies among different results: in our 
case, an inspection of the likelihood maps obtained separately 
on each field shows that the NTTDF, having two bright objects 
(Y ~ 25.5 - 25.7, approximately Mi^oo < -21.2 at z = 6.8), has 
a great effect in skewing the global likelihood towards brighter 
values of M, . We also note that som e theoretical models (e.g. 
iTrenti et alJl2010l:lFinlator et alJl20l"ol) predict a dimming of M» 
with redshift. However, several model parameters are largely un- 
constrained by the observations, while a large dust extinction 
might be r equired to match o bserved and predicted LFs at the 
bright end (lLacevetal.ll2010l) . 



6.4. Cosmic Variance 

The eff'ects of cosmic variance are reduced in our case, since 
our data come from three independent areas, albeit of dififer- 
ent sizes (the GOODS-South field being covered by two of 
the four Hawk-1 pointings). We evaluate the possible impact of 
cosmic variance using the mock catalogues of the Millennium 
Simulation ( Kitzbichler & White 2007) in the same way as dis- 
cussed in CIO. For each of the three Hawk-1 areas (GOODS- 
South, BDF, NTTDF), we extract 200 fields of the same size 
from independent Millennium light-cones, and we apply a cor- 
responding photometric selection criteria on galaxies at 6.5 < 
z < 7.4 (bracketing the peak of our selection window), with- 
out any constraint on the distribution of host haloes. We esti- 
mate that a cosmic variance of ~ 21% affects the total number 
counts of z-drop LBGs in our survey. We find that the evolution 
is still confirmed at a > 99% confidence level by our maximum 
likelihood approach even allowing a ~ 2 1 % variation in the to- 
tal number density. Indeed, after accounting for all the observa- 
tional effects, we estimate that we would have observed ~ 30 
z-drops in our survey in the case of a non-evolving LF: a factor 
of two higher than the observed number. However, while cosmic 
variance has not a significant effect on our conclusion that the 
LF strongly evolves from z ~ 7 to z ~ 6, it can have a great ef- 
fect in determining the form of this evolution. Cosmic variance 
is strongly luminosity dependent, and it is as high as 41% for 
galaxies brighter than Y = 25.8 in our survey, thus affecting the 
determination of the M, parameter 



Table 5. Properties of the z ~ 7 population" 



Puv 
SFRD 

log{Nu,n) 



1.5!;;^ 10-^ erg S-' Hz'' Mpc--' 



a - LFs in the 95% c.l. region (M1500 < -19.0) 



6.5. UV luminosity density, SFRD and constraints on cosmic 
reionization 

While the M, and (p parameters are highly degenerate, the num- 
ber density of bright galaxies, i.e. the integral of the bright end 
of the LF is much better constrained, and so are derived inte- 
gral quantities such as the UV luminosity density {puv) and star 
formation rate density (SFRD). 

We conservatively consider the model LFs within the 95% 
c.l. region of our likelihood analysis to derive the puv by inte- 
grating L ■ 0(L) up to the luminosity corresponding to Mi 500 = 
-19.0. We conve rt these value s in a SFRD following the stan- 
dard formula bv iMadau et al. (fl998) and applying the extinc- 
tion correction of Meu rer et al.l {1999) (considering an average 
UV slope /3 = -2.0). Finally, we use puv to evaluate the emis- 



sio n rate Njnn of hydrogen ion izing photons per Mpc^ follow- 



ing IBolton^&Haehnelg ^20071). We consider an escape fraction 
fesc = 0.2, a spectral index a, - 3.0 and an ionizing emission 
density at the Lyman limit eg = puyl6.Q. 

We report in Tab. |5]the range of values for puv, SFRD and 
log{Nio„). These values are perfectly consistent with the anal- 
ogous ones presented in CIO and derived from the LFs in the 
68% c.l. region of the GOODS li kelihood. Con s iderin g the same 
integral of the z=6 UV LF of McL ure et aP ([2009), our esti- 
mated Puv implies a drop of a factor ~ 3.5 in the UV lumi- 
nosity density from z=6. The lower limit for the ionization rate 
requir ed to balance recombination at z = 7, computed accord- 
ing to lMadau et al.l ( 1999) and assuming an HII clumping factor 
equal to one, is log(Nrec) - 50.1, which is a factor of two higher 
than the highest value allowed by our analysis. This demon- 
strates that, under usual assumptions, bright UV galaxies alone 
cannot keep the universe reionized at z ~ 7. By varying the es- 
cape fraction, we obtain that values larger than fesc = 0.5 are 
required to reconcile the emission rate from bright galaxies with 
the one required for reionization. This demonstrates that, either 
UV bright galaxies at z ~ 7 have different physical properties 
with respect to lower redshift LBGs, or, most probably, a crucial 
contribution to the reionization process comes from galaxies at 
the faint end of the LFs or from other kind of sources. Once dif- 
ferent integration limits are taken into ac count, our es timates are 
in agreement with the results obtained by Bouwens et al. (2008!) ; 
lOuchi et al.1 (l2009al) : iGonzalez et al.l (I201oIk 



7. Constraints on the LBG number density at z ~ 8 

We exploited deep Hawk-I J- and A'-band observations to put an 
upper limit on the number density z > 7.5, F-drop galaxies in our 
survey. We used the observations of the BDF and NTTDF fields 
presented in Sect. |2] and deep observations of the two GOODS- 
South pointings obtained both in our program as well as through 
a similar ESO observing program (Clement et al. in preparation). 
We obtained a multicolour catalogue with the 7-band as detec- 
tion image using SExtractor in dual mode over the full imag- 
ing set presented in this paper and in CIO. We used the same 
detection parameters and 2FWHM apertures adopted for the Y- 
detected catalogue, computing aperture corrected total magni- 
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Table 6. LBG number density at z ~ 8 

Mag. Range fiW^IpF^ 
Misoo < -22.0 <0m. 



tudes through appropriate corrections in each band. We chose 
the colour selection criteria in order to isolate galaxies having 
the Lyman-break sampled by the Y - J colour, and to exclude 
contamination from lower redshift galaxies on the basis of the 
expected colours for passive and dusty-starburst galaxies mod- 
elled as described in Sect. 13.21 



(Y-J) > 0.8 

(Y-J) > l.l+0.6-(J-K) 



We also required no detection in the optical bands adopting 
the same S /N criteria outlined for the selection of z-drop galax- 
ies. We limited our selection to 7=24.5,24.8,25.0 in the BDF, 
NTTDF and GOODS pointings respectively, up to which we es- 
timate that our catalogues are 100% complete, and we used the 
same area chosen for the selection of z-drop galaxies in order to 
avoid the noisiest regions in any image. 

With these criteria we found no candidate F-drop galaxy in 
our survey. Considered the average J - M1500 relation at the me- 
dian redshift of our colour selection (z=8), we are probing the 

M < M, region of the LF at M1500 22.5. We report in Tab. |6] 

an upper limit on the number density of very bright F-drop LBGs 
estimated as the inverse of the volume sampled by our survey in 
the redshift interval 7.5 < z < 9.0. 



8. Summary and conclusions 

We presented in this work the results of a F-band survey of the 
two high galactic latitude BDF and NTTDF fields aimed at de- 
tecting galaxies at z > 6.5 and measuring their number density. 
The survey is based on deep observations obtained under a ded- 
icated ESO Large Programme. We made use of Y, J, K band 
observations performed with Hawk-I, the new near-IR camera 
installed at the VLT, and of F0RS2 Z-band observations. We 
matched and combined these data with deep archive FORS 1 and 
FORS2 observations in the U, B, V, R, I filters to detect high 
redshift LBGs under the main criterion Z -Y > 1, requiring no 
optical detection and flat Y-J and Y - K colours. The colour 
selection criteria have been tailored in order to exclude lower 
redshift passive galaxies and dusty starbursts. Galactic T-dwarfs 
and galaxies exhibiting large Z - Y colours as well as signifi- 
cant emission in the optical bands, possibly intermediate redshift 
sources with bright emission lines. 

As a result, we isolated 8 highly reliable z-drop candidates 
in the magnitude range Y ^ 25.5 - 26.5 over a total area of 
70.1 arcmin^. We combined this z-drop sample with the similar 
one extracted from two pointings over the GOODS-South field 
comprising seven galaxies atY < 26.1 . 

We performed a stacking analysis of the 15 objects to es- 
timate the average properties of M ~ M» galaxies at z > 6.5. 
The photometric redshift of the stacked object is z = 6.85^^ j^ 
in perfect agreement with the estimated selection window of our 
survey. The stacked SED is fitted by an E{B - V) = 0.05 ^J},',^ at 
a 68% confidence level, indicating a low dust content in agree- 
ment with previous analysis of z~ 7-8 objects ^ Finkelste in et al.l 
l2009HSchaerer & de Barrosll2010HGonzalez et al.„20ia) . 



We then estimated the number density and the LF evolu- 
tion on the basis of detailed MonteCarlo imaging simulations 
accounting for all the uncertainties involved in the observations: 
detection completeness, photometric scatter, and random fluctu- 
ations in the S/N measure due to overlapping unresolved sources, 
or other eff'ects. We first computed a binned estimate of the 
galaxy number density at z ~ 7 following two different proce- 
dures. The first one, which is based on an average F-Misoo rela- 
tion and on an estimate of the redshift dependent completeness 
of our selection, is the procedure commonly adopted in the lit- 
erature. The second method is more conservative, and exploits 
a x^ minimization to compare the observed number counts to 
those predicted on the basis of MonteCarlo simulations for dif- 
ferent combinations of galaxy densities. This second procedure 
intrinsically considers the uncertainties in the Y-Mx^qq conver- 
sion due to photometric scatter, to the redshift distribution and to 
the intrinsic properties of different galaxy models. We find that 
the two procedures are consistent and they are in agreement with 
similar analysis from the literature. However, the more conser- 
vative procedure highlights that sources of statistical uncertainty 
are usually underestimated. 

To assess the degree of evolution of the UV LF at z > 6.0, 
we also simulated galaxy populations following different UV 
Schechter functions with linearly evolving parameters log{(p) 
and M». For each of the four Hawk-I pointings we compared 
the resulting distributions of simulated magnitudes and colours 
with the observed ones following a maximum likelihood ap- 
proach. We find strong evidence of evolution of the LF above 
z=6: our analysis rules out at a > 99% confidence level that the 
LF remains constant in both (p and M, above z = 6. Our likeli- 
hood maps for the Schechter parameters indicate a strong evo- 
lution in and a brightening of M, with redshift. However, the 
detection of two bright objects {Y ~ 25.5 - 25.7, correspond- 
ing to M1500 ^ -21.2) in the NTTDF pointing have a major 
role in skewing the evolution of M, towards bright values. The 
two Schechter parameters are, however, highly degenerate and 
our findings are also consistent within the uncertainties with a 
milder evolution of (p and a constant or s lightly dimmi ng M, 
as indi cated by other authors dBouwens et al...2008; .Ouchi et al.l 
l2009ah . We estimate that the possible eff'ect of cosmic variance 
is not capable of reconciling the observed number density of 
z-drop galaxies with the one predicted for a non-evolving LF. 
However, the strong dependence on luminosity of the cosmic 
variance, and the relatively small magnitude range probed by our 
survey at M < M», can influence the determination of the form 
of the evolving LF and provide an explanation for the difference 
between the evolution we determine and other estimates in the 
literature. 

The uncertainty and the degeneracy in the M* and best-fit 
values are not reflected in a comparable uncertainty in the num- 
ber density of bright galaxies. We conservatively consider the 
model LFs within the 95% c.l. region of our likelihood analysis 
to derive for galaxies at M1500 < -19.0 an UV luminosity den- 

' //z"i Mpc-' 



sity puv 



1. 5^1110^^ erg s-^ 



-0.8 ' 



a star formation rate 



density SFRD 



o 9+3.6 
-'■'^-1.9 



10 Mq yr Mpc ^ and an emission 



rate of hydrogen ionizing photons log{Nion) - 49.4;^q3 Mpc^^. 
The UV luminosity density is lower than the corresponding one 
at z ~ 6 by a factor ~ 3.5, while Nion is lower by at least a 
factor of ~ 2 than the l o wer U mit required for reionization ac- 
cording to 'Madau et al.l (Il999l) . considering f^sc - 0.2 and an 
Hll clumping factor equal to one. This implies that UV bright 
galaxies alone cannot reionize the universe, unless their physi- 
cal parameters are much different from those of lower redshift 
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LBGs (e.g. fcsc > 0.5, harder UV spectrum etc.). Most probably, 
the crucial contribution to reionization comes from galaxies at 
the faint end of the LF or from other kind of sources. Finally, we 
exploit the Hawk-1 J and K band observations of our survey to 
derive an upper limit of 2 ■ IQ^^Mpc^^ for the number density 

of M 22.5 LBGs at z ~ 8 from the non-detection of F-drop 

galaxies up to 7 ~ 25. 
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